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ABSTRACT 

The most massive Galactic globular cluster uj Cen appears to have two, or perhaps 
more, distinct main sequences. Its bluest main sequence is at the centre of debate be- 
cause it has been suggested to have an extremely high helium abundance oiY'^ 0.4. 
The same helium abundance is claimed to explain the presence of extreme horizontal 
branch stars of tu Cen as well. This demands a relative helium to metal enrichment of 
AY/AZ ~ 70; that is, more than one order of magnitude larger than the generally ac- 
cepted value. Candidate solutions, namely, AGB stars, massive stars, and supernovae, 
have been suggested; but in this study, we show that none of them is a viable channel, 
in terms of reproducing the high value of AY/AZ for the constrained age difference 
between the red and blue populations. Essentially no populations with an ordinary 
initial mass function, including those candidates, can produce such a high AY/AZ 
because they all produce metals as well as helium. As an alternative, we investigate 
the possibility of the stochastic "first star" contamination to the gas from which the 
younger generation of uj Cen formed. This requires the assumption that Population 
III star formation episode overlaps with that of Population II. While the required 
condition appears extreme, very massive objects in the first star generation provide a 
solution that is at least as plausible as any other suggestions made before. 
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I INTRODUCTION 

The most massive Galactic globular cluster, lj Cen, shows 
a double (or more) main sequence in the colour-magnitude 
diagrams (CMD), with a minority population of the bluer 
main sequence (bMS) and a majority population of red- 
der MS {rMS) stars (Anderson 1997; Bedin et al. 2004). 
The interpretation of the bA4S by a huge excess in the he- 
lium abundance {Y ~ 0.4) has been suggested via CMD 
fittings (Norris 2004) and supported by the spectroscopic 
study on the bMS stars (Piotto et al. 2005). The extreme 
helium abundance is claimed to explain the extremely-hot 
horizontal-branch (HB) stars of uj Cen as well (Lee et al. 
2005). 

The helium-enrichment parameter AY/AZ required to 
explain the sub-population of uj Cen are ~ 70, more than 
an order of magnitude larger than the currently-accepted 
value, 1-5 (Fernandes et al. 1996; Pagel & Portinari 1998; 
Jimenez et al. 2003). To explain this, several possible pol- 
luters, i.e., AGB stars, stellar winds associated with mas- 
sive stars during their early evolutionary phases and type 

II supernovae (SNe II) have been suggested (e.g. Norris 
2004; DAntona et al 2005). However, Bekki & Norris (2006) 
show that these candidates cannot reproduce such a high he- 



hum abundance for reasonable initial mass functions (IMFs) 
within the scheme of a closed-box self-enrichment. 



The acute point is not the high helium abundance itself 
but the high value of AY/AZ, because populations con- 
taining those candidates produce a large amount of metals 
as well as helium according to canonical stellar evolution 
models. Maeder & Meynet (2006) have recently suggested 
that "moderately" fast-rotating massive (~ 60 Mq) stars 
expel helium-dominant gas into space and thus could pro- 
vide a solution to this problem. We show that this scenario 
provides a viable solution only if additional conditions on 
the mixing and escaping of metals are invoked. 



As an alternative solution, we investigate the possibility 
that the late generation of "the first stars" (Population III) 
might affect the chemical mixture of the minority of Popu- 
lation II stars. We are particularly inspired by the work of 
Bromm & Loeb (2006) that suggested an overlap between 
the Population III and II star formation episodes. The chem- 
ical yields of such supposedly-heavy, zero-metal stars have 
been computed by Marigo, Chiosi, & Kudritzki (2003). We 
present in this Letter the result of our investigation. 
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2 CHEMICAL EVOLUTION MODELS 

In order to investigate the pattern of the hehum enhance- 
ment in a population, we employ a simple chemical evolution 
code and a realistic initial mass function (IMF). We describe 
a simple chemical enrichment model essentially following the 
formalism of Tinsley (1980) and reduce a set of a few pa- 
rameters following the formalism of Forreras & Silk (2000, 
2001). A two-component system is considered, consisting of 
cold gas and stellar mass. The net metallicity Z and helium 
content Y of two systems are traced. We assume instanta- 
neous mixing of the gas cjocta from stars and instantaneous 
coohng of the hot gas component. The mass in stars, Ms{t) 
and in cold gas, Mg{t), are normalised to the initial gas 
mass, 



M,{t) 



(1) 



Mg(0)' - M,{0) 

where the initial states of the galaxy is assumed to be com- 
pletely gaseous without stars: /tts(O) = 0. 

A Schmidt-type star formation law (Schmidt 1963) is 
assumed: 



v(t) = affM;(t) 



(2) 



where the parameter CcB implies the star formation effi- 
ciency. We assume n = 1, that is, a linear law ^{t) = 

CeffMg(t). 

Exponential infall of primordial gas is assumed: 

.p-{t-nag)/n„/ 



f{t) = e(t - riag)A,,,fe- 



(3) 



where, Q{t) is a step function. The parameters A^nf: ^inf^ 
and Tiag are the infall rate, timcscalc, and delay, respectively. 
In order to explain the G-dwarf problem (Tinsley 1980), gas 
infall has been considered eis a solution (Larson 1972). 

Some cold gas is heated to high temperature by su- 
pernovae and/or AGN and can be driven out: outflows. It 
is also an important factor to the final chemical properties 
of galaxies (Larson 1974; Arimoto & Yoshii 1987). We use 
a free parameter Bo-ut that represents the fraction of gas 
ejected following the formalism of Fcrrcras & Silk (2000). 
This parameter should be a function of the mass of the 
galaxy, whose potential well determines whether the winds 
are strong enough to escape its gravitational potential. 

There arc five input parameters : Ccff , -Bout, ^in/; '^inf-, 
Tiag with the following initial conditions: the initial metal- 
licity Zo = 10"'' and the initial helium of stars and gas 
Yo = 0.235, the IMF slopes and cutoffs. 



2.1 Mass Evolution 

The evolution of gas mass is given by 



^ = (1 - Bout)E{i) - C,slJ-g{t) 

+e{t - riag)Air,fe-^*-^^"9)/nnf 

where the stellar mass ejecta E{t) is defined as 

E{t) = / dm.(f>{m){m - Wm)CefllJ,g{t - Tm)- 

J mt 



(4) 



(5) 



where 0(m) is the Scalo IMF (Scalo 1986) with cutoffs at 
0.1 and IOOMq. We adopt the broken-power law of Ferreras 
& Silk (2000, Eq. 6) for the stellar lifetime. Remnant mass 



(6) 



(7) 



for a star with main sequence mass m, Wm, is adopted from 
Ferreras & Silk (2000). 

2.2 Chemical Evolution 

The chemical evolution of gas is given by 

= -C,aZg(t)tig(t) + {1- Bout)Ez{t) 

+@{t - riag)ZiAi^ie-(*-^^-9)/-ir.f 

Ez{t) = dm4>{m)CeB [{m - Wm){Zg^,g){t - Tm) 

+mpml^g{t - Tm)] , 

where pm denotes the mass fraction of a star of mass m 
that is newly converted to metals or helium and ejected. 
We approximate it by a polynomial fit to the pm prediction 
of Maedcr (1992). In Figure 1, we show Macdor's chemical 
yields for metal-poor {Z = 0.001) stars and our fitting func- 
tions. For our reference model, we use these chemical yields. 

We also show the metal yields of extremely metal-poor 
{Z = 0.00001) rotating stars (Maeder & Meynet 2006) 
in dotted line. Maeder & Meynet computed the yields for 
60 Mq stars, and in order to domonstrate their effects to 
the AY/AZ of an integrated population, we set up an ex- 
treme assumption that the newly-proposed metal yields may 
be applicable to all heavy stars of M ^ 50. 

The crucial point shown here is that canonical stellar 
models produce and spread into space metals as well as he- 
lium. No combination of these stars can produce AY/AZ 
that exceeds the values of the constituent stars. The rotat- 
ing metal-poor massive star models by Maedcr & Meynet 
(2006) arc exception and indeed can produce high values of 
AY/AZ. So we explore below whether a population includ- 
ing such stars can indeed present a solution to our problem. 

SNe la are generally considered important in chemical 
evolution studies because they eject a considerable amount 
of iron-peak elements: ~ 0.7 Mq (Tsujimoto et al. 1995). 
For our exercise, we have tried the most up-to-date SN la 
rate of Scannapieco & Bildsten (2005), consisting of two 
components: star formation rate and the total stellar mass 
of the system. The contribution by SNe la is negative to 
our study because they are copious producers of metals. We 
remove its effects for simplicity. 

We have found that gas infall or outflows of the pro- 
cessed material have negligible impacts on AY/AZ unless 
specific chemical elements alone (e.g., helium or metals) are 
affected. Hence, only the closed-box model is considered. 

The star formation efficiency parameter CcH has a small 
effect. When a range CcH = 0.1 - 5.0 Gyr^^ is considered, 
a population can reach AY/AZ ~ 1 - 2 after a 1 Gyr evo- 
lution. We show a typical case of Cgs = 1 in Figure 2 as 
a solid line. The value of our reference model is consistent 
with the most recent observational mean value of the helium 
enrichment parameter of the galactic gas. 

Having failed to recover the high AY/AZ in the stan- 
dard chemical enrichment model, we consider the stochastic 
effect in the helium yields pm,Y. We assume that the stars 
in a specific mass range can have larger values of Pm,Y than 
the current stellar physics suggest. 

Fig. 2 presents part of the result of the stochastic-effect 
test, where Cgg = 1 is assumed; (1) when Pm,Y is doubled 
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Figure 1. top: The chemical yields of metal-poor (Z = 0.001) 
stars from Maeder (1992) and our fits to them. The red dash-dot 
line is for our ad-hoc case employing the Z = 0.00001 rotating 
massive models by Maeder & Meynet (2006) for M ^ SOMq. 
bottom: The resulting helium-to-metal ratios (bottom). The red 
line is for the "ad- hoc case" . 



< 



>- 
< 



Normol system 
doubled p„ 
AGB 
SN II 

massive rotating stars 




2 4 6 8 10 12 14 
t/Gyr 

Figure 2. The helium enrichment parameters resulting from the 
stochastic effect. Our reference model is shown in solid line. Other 
lines are for the cases that p„ y is doubled for (1) all stars, (2) 6 
- 10 Mq stars (AGB progenitors) stars, and (3) 40-50 Mq (SN II 
progenitors). Finally, triple-dot dashed line shows the model that 
assumes that all heavy stars (M ^ 50) are fast-rotating showing 
yields of Maeder & Meynet (2006). All models are for the case of 
Cefi = 1 and the chemical mixing time scale of 50 Myr. 



for all the stars, (2) when Pm,Y is doubled for 6-10 Mq 
stars, maximising the AGB effect, and (3) when Pm,Y is 
doubled for 40 - 50 Mq stars, allowing the SN II effect to 
dominate. Lastly, the triple-dot dashed line is for the model 
where all heavy stars (M ^ 50) are fast-rotating showing 
yields of Maeder & Meynet (2006). Because instant recy- 
cling of chemical elements is unrealistic, we use the chemical 
mixing time scale of 50 Myr. 

The Ay/AZ after a 13 Gyr evolution of the AGB- 
dominant system is merely ~ 2 times greater than that of 
the reference system. The effect is even smaller for the SNe 
II dominant system. Even if we double the values of Pm,Y 
for all stars, the system could attain only l\Y / l^Z « 4. But 
of course, for our purpose, we need a much higher helium 
enrichment within a much shorter timescale (« 1 Gyr). 

Doubling the helium yields is ad hoc and hard to justify 
in terms of nucleosynthesis. Even in this unphysical test, 
however, the stochastic effect does not help us achieve the 
high Ay/AZ. This is because (1) ordinary stars produce 
and spread into space not only helium but also metals and 
(2) most notably in the case of the "Maeder & Meynet case" , 
there are many more low-mass stars in the ordinary IMF 
than massive stars. 

Other possible contributing factors, e.g., IMF slope and 
the remnant mass Wm, also have some influence. One may 
wonder about the effect of radically-different IMFs (e.g., 
with a reversed slope). This indeed would produce more he- 
lium but even more metals in ordinary stellar populations, 
hence in fact lowering the resulting value of ISY //S.Z. On 
the other hand, a reversed IMF slope would dramatically in- 
crease Ay/AZ in the "Maeder & Meynet case". However, 
this should not be common because otherwise all popula- 
tions with massive stars should have enhanced helium. 



When assuming that the helium enrichment of hMS of 
UJ Cen is due to the ejecta from the rMS formed earlier, 
the candidate solution is even more severely constrained. 
For example, low-mass stars which produce mainly helium 
can generate high values of Ay/AZ upto ~ 8 but they take 
much too long to spread the processed materials into space, 
when compared to the age difference between the two sub- 
systems in UJ Cen, i.e., 1-2 Gyr (Lee et al. 2005). Any star 
whose lifetime is longer than this would have no impact at 
all. 

The "Maeder & Meynet case" is noteworthy as it is the 
only scenario where high values of Ay/AZ can be achieved 
during a short period of time. Assuming all stars form si- 
multaneously in a population, Ay/AZ can be as high as we 
need during the first few Myr until intermediate-mass stars 
begin to spread a large amount of gas into space. In this 
sense, this scenario provides a solution if the bMS popula- 
tion of UJ Cen is only slightly younger than its rMS popu- 
lation. However, the time window {^^ 10^~* yr) for this to 
happen is very small compared to the constraint from CMD 
studies (Lee et al. 2005). The possibility of difference in the 
effect of stellar wind for helium and metals demands further 
studies, and the role of the local black hole accretion on 
metals may also be noteworthy (Maeder & Meynet 2006). 



3 VERY MASSIVE OBJECTS 

We finally explore the plausibility of the "first" Pop HI 
star pollution. The first stars in the universe are thought 
to have formed out of the metal- free gas at redshifts z>10 
(Bromm & Larson 2004; Ciardi & Ferrara 2005). These Pop 
III stars are often predicted to be very massive, M > IOOA^q 
(Bromm, Coppi & Larson, 2002; Abel, Bryan & Norman 
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Figure 3. (Top) Filled circles represent the Pm.i'i Vm,Zi newly 
synthesized and ejected helium or metal mass fraction, calculated 
from yields of Marigo et al. (2003) for 120, 250, 500, 750 and 1000 
Mq VMOs. Dotted line implies linear interpolation of given 5 
points in logarithmic scale. (Bottom) The AY/AZ of each mass 
of VMOs 



2002), although a different suggestion has also been made 
(Silk & Langer 2006). They are considered to be significant 
contributors to various phenomena (Carr et al. 1984). 

We have calculated the chemical enrichment evolution 
of the closed-box system made only of VMOs. Marigo et al. 
(2003) give the evolutionary properties of VMOs including 
lifetime, helium and metal yields, carbon and helium cores at 
central C-ignition and mass loss for I2O-IOOOM0 . Two dif- 
ferent mass loss rates are considered by them: the radiation- 
driven mass-loss model and the rotation-driven model. The 
yields and mass loss rates of the two models are slightly 
different from each other, and we use the rotation-driven 
model because it produces a larger amount of helium. We 
find linear fits to these yields in the logarithmic scale (Fig. 3) 
and show the resulting AY/AZ for each mass of VMO. The 
helium enrichment parameter of VMOs ranges 63 - 6 x 10^! 

We have calculated AY / /S.Z for a population of VMOs 
with various mass cuts and IMF slopes. The l^Y / AZ value 
after the 1 Gyr evolution in a closed-box system varies with 
IMF slope and mass cut. The values of AF/AZ as a func- 
tion of lower mass cut are shown in Fig 4. We fix the high 
mass cut as lOOOAfQ assuming the Salpeter IMF. The re- 
sulting range of AF/AZ after IGyr of evolution is two 
orders of magnitude higher than that of an ordinary sys- 
tem. Our model matches the values of both Y (~ 0.4) and 
Ay/AZ (~ 70) at age ~ 1 Gyr suggested for the hMS sub- 
population, as marked in Fig. 5. Interestingly the solution 
suggests M{VMO) ~ 995-1OOOM0 ; but considering the un- 
certainty in the VMO yields, it may not be significant. 
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Figure 4. The mass-cut dependence of AY/AZ of the ejecta 
from VMOs. We constrain high mass cut of IMF as 1000 Mq and 
IMF slope as canonical value, x = 1.35. The system consisting of 
12O~lOOOA/0 VMOs with Salpeter IMF slope results in AY/AZ~ 
520 and 995~1OOOM0 system produces AY/AZ~ 72. 
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Figure 5. Chemical evolution in the gas of the VMO system 
for different values of star formatin efficiency. The values of he- 
lium and AY/AZ of the bMS population of u) Cen are marked 
(squares) in each panel. 



4 DISCUSSION 

By means of a simple chemical evolution test, we conclude 
that ordinary populations cannot produce AY/AZ much 
greater than 4. None of the candidates succeeds in repro- 
ducing the large helium enrichment claimed for the blue MS 
stars of uj Cen. The reason for this is simple: the evolution of 
ordinary stars result only in a modest value AY/AZ ~ 1-5. 
We confirm the result of Karakas et al. (2006) who rejected 
the AGB solution on a similar ground. The massive rotating 
star models of Maeder & Meynet (2006) reach high values of 
AY/AZ, but it is difficult to maintain the high value longer 
than a few Myr in a population with a realistic mass spread. 
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We present Pop III very massive objects (VMOs) as an 
alternative solution. They produce both the helium abun- 
dance and the helium-to-metal ratio searched for at the 
age t Ri 1 Gyr demanded by the observational constraint. 
This scenario requires an overlap between the star formation 
episodes of Pop III and Pop II. Suppose that lj Cen formed 
in the era of overlap of Pop III and Pop II formation. Ear- 
lier generations of Pop III stars form and their mass ejecta 
are mixed in the short timescale of roughly 10^~* yrs. First 
generation of Pop II stars start forming globally (this is not 
part of u> Cen) and the processed gas gets recycled and mixed 
in the proto-system until the gas reaches <Z>^ 0.001. By 
now, the mixed mean helium abundance is not any more 
extraordinary but with severe irregularities. The rMS pop- 
ulation of u) Cen forms with the mean chemical composi- 
tion. A group of Pop III VMOs form in the tail of the Pop 
III episode nearby out of a pristine (irregularity) gas cloud. 
Soon after this, that is, before their ejecta are mixed fully, 
the younger stars of w Cen (bMS) form under the heavy 
influence of this helium-rich ejecta. After the dynamical re- 
laxation, the two chemically-distinct populations could be 
in one system. Since this is a geographic and chronological 
stochastic effect, it would not be significant in the galactic 
scale. 

The plausibility of our scenario strongly depends on the 
validity of Marigo et al.'s yields. Other groups have com- 
puted the yields for zero-metallicity VMOs (e.g. Bond et al. 
1984; Ober et al. 1983; Klapp 1984), and Marigo et al.'s 
yields are not in a smooth continuation with the yields of 
low-mass zero-metal stars given by Limongi & Chieffi (2002) . 
The discontinuity may not be a big issue if the mass loss in 
VMOs may happen in an extreme fasion (Smith 2006). Con- 
sidering the acute dependence of our conclusion on the relia- 
bility of the chemical yields, it is urgent to perform detailed 
independent stellar evolution modelling for VMOs. 

The second question about this scenario is on the con- 
dition of VMO formation. The Pop II star formation is al- 
lowed in the broad regions of Pop III objects (Tsujimoto et 
al. 1999; Susa & Umemura 2006). The bMS population of 
ui Cen has ^ 30% of the cluster mass (Lee et al. 2005) . As- 
suming Mtot ~ 10® Mq, the helium-rich population mass 
would be MbMS ~ 3 x 10^ Mq. Assuming the primor- 
dial helium abundance is 23%, 17% of the bMS population 
mass, Mhb ~ 5 x 10** Mq, is the newly-generated helium 
mass. To generate 5 x 10** Mq pure helium from 1000 Mq 
VMOs, we need at least ~ 200 VMOs simultaneously. How- 
ever, Abel, Bryan, & Norman (2002) suggest that metal- 
free stars form in isolation due to the immense radiation 
from the first-forming star. Moreover, such a large number 
of Pop III stars forming close to each other, if that is pos- 
sible at all, would make the site extremely hostile for the 
younger generation of uj Cen to form. This poses serious 
challenges to our scenario. By definition, this stochastic ef- 
fect works only on small scales and not in the galactic scale. 
But the number of VMOs required would be smallest in the 
deepest local gravitational potential, which means that this 
rare event of achieving a high AF/AZ, would prefer larger 
gravitationally-bound objects, such as massive globular clus- 
ters rather than small ones. 
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